CONTROL OF UNDERWATER VEHICLES IN INVISCID FLUIDS. 

I: IRROTATIONAL FLOWS. 



RODRIGO LECAROS^'2 AND LIONEL ROSIER^ 



Abstract. In this paper, we investigate the controllability of an underwater vehicle immersed 
in an infinite volume of an inviscid fluid whose flow is assumed to be irrotational. Taking as 
control input the flow of the fluid through a part of the boundary of the rigid body, we obtain a 
flnite-dimensional system similar to Kirchhoff laws in which the control input appears through 
both linear terms (with time derivative) and bilinear terms. Applying Coron's return method, 
we establish some local controllability results for the position and velocities of the underwater 
vehicle. Examples with six, four, or only three controls inputs are given for a vehicule with an 
ellipsoidal shape. 

o 

c3 1. Introduction 

^ The control of boats or submarines has attracted the attention of the mathematical community 
|> from a long time (see e.g. [2l[3lll[9l[l0l[l6l[I71[l8l|19].) In most of the papers devoted to that 
issue, the fluid is assumed to be inviscid, incompressible and irrotational, and the rigid body 
(the vehicle) is supposed to have an elliptic shape. On the other hand, to simplify the model, 
the control is often assumed to appear in a linear way in a finite-dimensional system describing 
^ the dynamics of the rigid body, the so-called Kirchhoff laws. 

-t-j A large vessel (e.g. a cargo ship) presents often one tunnel thruster built into the bow to 
make docking easier. Some accurate model of a boat without rudder controlled by two propellers, 
, the one displayed in a transversal bowthruster at the bow of the ship, the other one placed at 
the stern of the boat, was derived and investigated in [T2|. A local controllability result for 
the position and velocity (six coordinates) of a boat surrounded by an inviscid (not necessarily 
irrotational) fluid was derived in |12j with only two controls inputs. 

The aim of this paper is to provide some accurate model of a neutrally buoyant underwater 
r vehicle immersed in an infinite volume of ideal fluid, without rudder, and actuated by a few 
number of propellers located into some tunnels inside the rigid body, and to give a rigorous 
analysis of the control properties of such a system. We aim to control both the position, the 
attitude, and the (linear and angular) velocities of the vehicle by taking as control input the 
^—i flow of the fluid through a part of the boundary of the rigid body. The inviscid incompressible 
sjjl fluid is assumed here to have an irrotational (hence potential) flow, for the sake of simplicity. 
• 'H The case of a fluid with vorticity will be considered elsewhere. 

rS Our fluid-structure interaction problem can be described as follow. The underwater vehicle, 
d represented by a rigid body occupying a connected compact set S(t) C M^, is surrounded by 
an homogeneous incompressible perfect fluid filling the open set Q{t) := \ §{t) (as e.g. for 
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a submarine immersed in an ocean). We assume that ^{t) is C°° smooth and connected. Let 
S = S(0) and 

n = 17(0) = R''^ \ S(0) 

denote the initial configuration (t = 0). Then, the dynamics of the fluid-structure system are 
governed by the following system of PDE's 

— + {u-V)u + Vp = (}, te{0,T), xen{t), (1.1) 

divn = 0, t e {0,T), X en{t), (1.2) 
u-u = {h' + ojx{x-h))-u + w{t,x), t e {0,T), X e dn{t), (1.3) 
lim u{t,x) = 0, te{0,T), (1.4) 

\x\—>-+oo 



rriQh" = J pvda, te{0,T), (1.5) 
dn{t) 

-{QJoQ*u)= J {x-h)xpvda, te{0,T), (1.6) 



dn{t) 

A. 

dn{t) 

Q' = S{uj)Q, tG(0,T), (1.7) 

u{0, x) = uo{x), X G Q, (1-8) 

ih{0),Q{0),h'{0),uj{0)) = {ho,Qo,hi,ojo) G x S0(3) x x (1.9) 

In the above equations, u (resp. p) is the velocity field (resp. the pressure) of the fluid, h 
denotes the position of the center of mass of the solid, u denotes the angular velocity and Q the 
3 dimensional rotation matrix giving the orientation of the solid. The positive constant mo and 
the matrix Jq, which denote respectively the mass and the inertia matrix of the rigid body, are 
defined as 

mo = J p{x)dx, Jq = J p{x){\x\'^Id — xx*)dx, 
§ S 

where p(-) represents the density of the rigid body. Finally, v is the outward unit vector to 
dQ{t), X X y is the cross product between the vectors x and y, and S{y) is the skew-adjoint 
matrix such that S{y)x = y x x, i.e. 

s{y) = 

The neutral buoyancy condition reads 

= ,1.10, 

s s 

/' (or /) stands for the derivative of / respect to t, A* means the transpose of the matrix A, 
and Id denotes the identity matrix. Finally, the term w{t, x), which stands for the flow through 
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the boundary of the rigid body, is taken as control input. Its support wiU be strictly included 
in dQ{t), and actually only a finite dimensional control input will be considered here (see below 



(1.17) for the precise form of the control term w{t,x)). 

When no control is applied (i.e. w{t,x) = 0), then the existence and uniqueness of strong 
solutions to (1.1)-(1.9) was obtained first in |20j for a ball embedded in M?, and next in |21] 
for a rigid body S of arbitrary form (still in M?). The case of a ball in was investigated in 



and the case of a rigid body of arbitrary form in M'^ was studied in )25j . The detection of 
the rigid body §{t) from partial measurements of the fiuid velocity has been tackled in [5] when 
n{t) = no\ (^0 C being a bounded cavity) and in [6] when n{t) = \ S(t). 



Here, we are interested in the control properties of (1.1)-(1.9). The controllability of Euler 
equations has been established in 2D (resp. in 3D) in [7J (resp. in [H]). Note, however, that 
there is no hope here to control the motion of both the fluid and the rigid body. Indeed, 0(i) is 
an exterior domain, and the vorticity is transported by the flow with a flnite speed propagation, 
so that it is not affected (at any given time) far from the rigid body. Therefore, we will deal 
with the control of the motion of the rigid body only. As the state of the rigid body is described 
by a vector in M^'^, it is natural to consider a finite-dimensional control input. 

Note also that since the fluid is flowing through a part of the boundary of the rigid body. 



additional boundary conditions are needed to ensure the uniqueness of the solution of ( 1.1 )-( 1.9 ) 
(see [13], [H]). In dimension three, one can specify the tangent components of the vorticity 
C{t,x) := curl v{t,x) on the inflow section; that is, one can set 

C(t, x)-Ti = Co{t, x) ■ Ti for w{t, x) < 0, i = 1, 2, (1.11) 

where C,Q{t,x) is a given function and Tj, i = 1,2, are linearly independent vectors tangent to 
d^l{t). As we are concerned here with irrotational flows, we choose Co = 0. 

In order to write the equations of the fluid in a fixed frame, we perform a change of coordinates. 
We set 

X = Q{t)y + h{t), (1.12) 
v{t,y) = Q*{t)u{t,Q{t)y + h{t)), (1.13) 
q{t,y) = p{t,Q{t)y + h{t)), (1.14) 
l{t) = Q*{t)h'{t), (1.15) 
r{t) = Q*{t)u{t). (1.16) 

Then x (resp. y) represents the vector of coordinates of a point in a flxed frame (respectively 
in a frame linked to the rigid body) . We may without loss of generality assume that 

/i(0) = 0, Q(0) = Id. 

Note that, at any given time t, y ranges over the fixed domain when x ranges over il(t). 
Finally, we assume that the control takes the form 



m 

w{t,x) = w{t,Q{t)y + h{t)) = ^Wj{t)xj{y), 



(1.17) 
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where m G N* stands for the number of independent inputs, and 'Wj{t) G M is the control input 
associated with the function Xj £ C°°{dQ). To ensure the conservation of the mass of the fluid, 
we impose the relation 



J Xj{y)d(y = for 1 < j < m. 

an 



Then the functions (u, q, I, r) satisfy the following system 



(1.18) 



dv 
di 



+ {{v — I — r X y) ■ 'V)v + r x v + 'Vq = 0, 

div V = 0, 

v-v={l + rxy)-v+ ^ Wj{t)xj{y), 

l<j<m, 

lim v{t, y) = 0, 



qu da — niQr x Z, 



do. 



(/(0),r(0)) 



q{y X v)da — r X Jor, 
an 

-- (/ii,a;o), v{0,y) = uo{y). 



(0,T), y eO, 


(1.19) 


(0,T), yen, 


(1.20) 


(0,T), yedn 


(1.21) 


(0,T), 


(1.22) 


(o,r), 


(1.23) 


(o,r), 


(1.24) 




(1.25) 


p:Tl)-(1.9)by 


assuming that 



laws, in which the control input w appears through both linear terms (with time derivative) and 



bilinear terms. The investigation of the control properties of (2.65) is performed in Section 3. 
After noticing that the controllability of the linearized system at the origin requires six control 
inputs, we apply the return method due to Jean-Michel Coron to take advantage of the nonlinear 
terms in ( 2.65| ). (We refer the reader to for an exposition of that method for finite-dimensional 
systems and for PDE's.) We consider the linearization along a certain closed-loop trajectory 



and obtain a local controllability result (Theorem 3.11) assuming that two rank conditions are 
fulfilled, by using a variant of Silverman-Meadows test for the controllability of a time-varying 
linear system. Some examples using symmetry properties of the rigid body are given in Section 
4. 



2. Equations of the motion in the potential case 

In this section we derive the equations describing the motion of the rigid body subject to flow 
boundary control when the fluid is potential. 



2.1. Null vorticity. Let us denote by 

C{t,y) = cuvlv{t,y) := (V x v){t,y) 
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the vorticity of the fluid. Here, we assume that 

Co = curl vq = in 17 
and that the three components of C are null at the inflow part of di}, namely 

m 

C{t,y)=0, if y G Ui<j<r„ Supp Xj and ^ u;^ (t)xj(y) < 0. 



Proposition 2.1. Under the assumptions (2.1) and (2.2), one has 

C = cnilv = in [0, T] x 0, 



Proof. Let us introduce v := v — I — r x y. Then it follows from (1.20) that 

div(i>) = 0, 

and 

curl{v) = C — 2r. 
Applying the operator curl in ( |1.19 ) results in 



dt 



+ cml{{v ■ V)v) + curl(({; • V)(/ + r x y)) + curl(r x v) = 0. 



We note that the following identities hold: 

curl((t' • V)^) = {v ■ V)curl(t;) — (curl(t;) • V)^ + div(f )curl(w) 

and 

{v ■ V)(r X y) = r X V, curl(r x v) = div(w)r — (r • V)f . 
Using (2.4)-(2.8), we see that C satisfies 



dt 



+ {v-V)C-{C-V)v = 0. 



Let ip = ip{t, s, y) denote the flow associated with v, i.e. 

dip 



v{t,<p), with (p\t=s = y. 



2.1 
2.2 

2.3 

2.4 
2.5 

2.6 

2.7 
2.8 

2.9 
(2.10 



We denote by G{t,s,y) = j^{t,s,y) the Jacobi matrix of if. Differentiating in (2.10) with 



{t, ip{t, s, y)) • G{t, s, y), where G{s, s, y) = Id (identity matrix). 



(2.11 



respect to yj {j = 1, 2, 3), we see that G{t, s, y) satisfies the following equation: 
dG__dv^ 
dt dy 

We infer from ( p^ and \2.ll\ that 

det G(t,s,y) = 1. (2.12 

Following Yudovich [13], we introduce the time t*{t,y) G [0,t] at which the fluid element 
first appears in Q, and set y*{t,y) = ip{t*{t,y),t,y). Then either t* = 0, or t* > and y* S 
Ui<j<m suppxi C dn with Ejli^i(**)Xj(y*) < 0. Set f{s,t,y) = G'^{s,t,y)C{s,ip{s,t,y)). 
Prom (IZ9l)-(|2.12|), we obtain that 



ds 



{s,t,y) = 0. 



(2.13) 
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Finally, integrating with respect to s in (2.13) yields 

C{t,y) = G'\t*,t,y)C{t*,y*), 



(2.14) 



which, combined to (2.1) and (2.2), gives (2.3). The proof of Proposition 2.1 is complete. □ 



Remark 2.2. The issue whether the result in Proposition 2.1 still holds with (2.2) replaced by 

m 

C{t,y) • Ti = 0, z = 1,2, if y £ Ui<j<m Supp Xj and ^Wj{t)xj{y) < 0, 
seems challenging. We notice that the result in [14J was proved solely when 0, was a cylinder. 



2.2. Decomposition of the fluid velocity. It follows from (1.20), (1.22) and (2.3) that the 
flow is potential; that is, 

V = V$, 

where $ = ^{t, y) solves 

A$ = 0, in(0,r)xJ7, 
{l + ryiy)-v+ Wj{t)xj{y) on{0,T)xQ, 

l<j<m 



di 



(2.15) 

(2.16) 
(2.17) 



lim V$(t,y) = 0, on(0,r). 

|j/|-5>+00 



Actually, <I> may be decomposed as 

'^{t,y)= ^ {kepi + rupi} + ^ Wjipj 

l<j<3 l<i<m. 

where, for i = 1, 2, 3 and j = 1, m, 

A(pi = Aifi = Aipj = in $7, 

Xj on on, 



du 



lim V(j)i{y) = 0, lim Vipi{y) = 0, lim V^j(y) = 0. 

\y\^+co \y\^+oo |y|-s>+oo 



(2.18) 
(2.19) 

(2.20) 

(2.21) 
(2.22) 



As the open set Q and the functions Xji ^ ^ j ^ "t-j supporting the control are assumed to 
be smooth, we infer that the functions V<j)i {i = 1, 2, 3), the functions V(^i (i = 1, 2, 3) and the 
functions V'i/'j (1 < J < rn) belong to H°^{0,). 

2.3. Equations for the linear and angular velocities. For notational convenience, in what 
follows f^f (resp. fg^f) stands for f^f{y)dy (resp. fg^ f{y)da{y)). 

Let us introduce the matrices M,J,N £ M^xs^ C^^,C-^ G M^xm^ L^^ , L-^ , e M^^^ 



and the matrices , G 



p3xm 



for p G {1, m} defined by 



J, 



Vv9i • Vv3j = {y X u)i(pj 



Of 



dv 



(2.23) 
(2.24) 



an 
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Nij = J V4>i ■ Vifj = j ViLpj = j (pi{y xi^)j, (2.25) 
n an an 

(C^'kj = y V</>. • V^j = J ui^jj = j cPiXj, (2.26) 
n an an 

(C^kj = J"^^i- VV'i = ^(y X zy)iVj = j ^iXj, (2.27) 
n an an 



{L^\, = j{V<Pj)iXp, i^)i,, = J iyxV^,)iXp, (2.28) 
an an 

(Rpkj = j (Vy^.Oap, (Rpkj = jiyx '^Vj)iXp, (2.29) 
an an 

(w^kj = I (vVi)ai., {Wp'kj = jiy^ v^,)ap. (2.30) 
an an 
Note that M* = M and J* = J. 

Let us now reformulate the equations for the motion of the rigid body. We define the matrix 
8 G M*'"'^ by 

' mold \ / M N \ 

Jo J^\N* J )■ ^^-^^^ 

It is easy to see that J is a (symmetric) positive definite matrix. We associate to the (linear 
and angular) velocity (/, r) € x of the rigid body a momentum- like quantity, the so-called 
impulse (-P, H) G x M^, defined by 

3(') = (n)- 

We are now in a position to give the equations governing the dynamics of the impulse. 
Proposition 2.3. The dynamics of the system are governed by the following Kirchhoff equations 

^ + C''w = {P + C^^w)xr- WpiL'^H + R^'r + W^'w], 

l<p<m 

^ + C'^w = (u + c-^yj)xr + {P + C^^w)xl- Wp {1:1^1 + R^r + w} , 

l<p<m 

(2.33) 

where w{t) := {wi{t), ...,Wmit)) G M*" denotes the control input. 



Proof. We first express the pressure q in terms of l,r,v and their derivatives. Using (2.3), we 
easily obtain 



hp 



V ■ Vv = ^-iid (r X y) • Vv — r x v = V((r x y) ■ v) (2.34) 
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Thus ( |L19D gives 
-Vq = 



dv „ /bp , , , 
+ V — I ■ V — (r X y) ■ V 



V j ^ + nipi} + ^ + ^ I ■ v - {r x y) ■ v 

,l<j<3 l<j<m 



hence we can take 



l<i<3 



(/ + (r X2/))-t;l (2.35) 



l<jr<m 



Replacing q by its value in (1.23) yields 



^ i<i<3 ^ - - ^ y- y_ 



9f7 an i<i<m 



+ (r X y)) ■ vj V 
(2.36) 



Using (2.34) and (1.20)-(1.21), we obtain 



-V = / V 



an 



/ 

— y (divt;)t> + y (i 

j {{l + rxy)-iy)v+ ^ u;j(t) ^ Xj(y)^^- 



(V ■ I'jV 



an 



Using Lagrange's formula: 



we obtain that 



an 

Now we claim that 



a X (b X c) = {a ■ c)b — (a ■ b)c, Va, 6, c G M^, 

{{I + r X y) ■ u) V — {{I + r X y) ■ v) = j {I + r x y) x {v x v). 

an 

i/xV/ = 0, V/GC2(n). 



(2.37) 

(2.38) 
(2.39) 

(2.40) 



an 



To prove the claim, we introduce a smooth cutoff function pa such that 



pa{y) 



1 if \y\ < a, 
if \y\ > 2a. 
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Pick a radius a > such that S C B{0, a), and set 

f{y) = f{y)pa{y). (2.4i) 

Then 

V/(y) = V/(y), yyGdn, 
and using the divergence theorem, we obtain 

vxVf = j vxVf = j curl(V/) = 0. 

an dn n 



Therefore, using (|2.40|) with f = ^ where V<l> = v, we obtain 

I X (v X u) = 0. 



an 



Another appHcation of (2.40) with / = yields 



ViV XV = j ^ ^ ej$, 



(2.42) 



(2.43) 



where {61,62,63} denotes the canonical basis in M'^. It follows from (2.42), (2.43), and (2.38) 
that 

{I + r X y) X {v X u) = r X j 
do. dn 

= X / MyHy) + n{t) J ^.(yMy)] + f2wj{t) I V',(y)Ky))- (2-44) 

an an an 



Combining (2.36) with (2.37), (2.39), and (2.44) yelds 



3 /i ft fy\ ft 

an an an 

m I 3 „ „ m .. 

j=i [j=i p=i 

3 AT?/ ft 



an 



an 



an 



-mor X I. 



(2.45) 



Let us turn our attention to the dynamics of r. Substituting the expression of q given in 



(2.35) in (1.24) yields 
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Jq-T = —r X JqV — < 
l<j<3 



k J 4>i{y X z^) + f 
k an on 



J i^^J'^" an 



V'j(y X ^) 



an 



{l + {rxy))-v]{yx u) (2.46) 



Prom [15^ Proof of Lemma 2.7], we know that 

|?;(y)| = |V$(y)| = 0(|y|-2), \Vv{y)\ = 0{\y\-^) as |y| ^ oo, 



so that 



veL^{n), \y\-\v\-\Vv\e L^{n). 



Note that, by (2.34) and (1.20) 



(2.47) 



div( — (cixy)) = V{ — )-{eiXy) 

= (v ■ Vv) ■ {ei X y) 

= V ■ V(y X v)i 

= div{{y X v)iv), 



-div(ei X y) 



and hence, using (2.47) and the divergence theorem, 



an 



-{y X u) 



an 



-{ei xy)-u 



div (^(ei X y)) 



div [{y X v)iv) 
{v • u){y X v)i 



an 



j {l + rxy)-v{yxv)i+ ^ Wj{t) j Xjiv x v) 



(2.48) 



an 



l<j<m 



an 
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Furthermore, using (2.38) we have that 



{I + r X y) ■ u{y X v)i - (/ + (r x y)) ■ v{y x v)i 
{l + rxy)- ^((ci xy)-v)u- ((e^ x y) ■ u)% 

{l + rxy)- {{a xy)x{vx v)) . 



an 



Combining the following identity 

3 

y^(a X Cj) X (ej X b) = —{a x b), Va, b £ 



with (2.43), we obtain 

/ I ■ {{a xy)x{uxv)) = V / / • ((ei x ej) x {u x yjv)) 
JdQ ~^ J 



J=i an 

3 



^ / • ^(ci X Cj) X j {vx yjv)^ 
an 

(^{a X Cj) X j {ej X 1/)$^ 
■''=^ an 

-I- J {eiXiy)^ = J (Ix 



For any given / G C^(r2), let 



/ := J (rxy)- {{a x y) x {u x V/)) 
ao. 



f still denoting the function defined in (2.41), we have that 

3 » 

I = {rxy) - {{e^ x y) x {ej x V/)) 

i=i an 

E { /(^ X e,) • ((e. X y) x (e, x V/)) } + E { / 

i = l Q .7 = 1 n 



(2.49) 



(2.50) 



(2.51) 



(rxy)- (ci X ej) x {ej x V/) 



n 



+ E ^ / (r X y) • ( (e» X y) X (ej x djVf) 
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Using again (2.50), we obtain 

3 



Yj{ / (ei X y) • ( (^^ X Cj) X {ej X V/) 

+ j {r-xy)- [{ei X 2/) X rot(V/)^ 
{a xy)-{rx V/) - {rxy)-{eiX V/) 



(r X y) • (cj X V/) 



n 



r • ((ei xy)x V/) - y r • (y x (e^ x V/)) 
= - j r-^{eiXy)xVf + yx{eixVf)^ 

= - j r ■ \^ei X {y xVf)^ = j (^r x {y x ^f)). = j {r x {y x v)f)^ , 
n n an 

where we used Jacobi identity 

a X {b X c) + b X {c X a) + c X {a X b) = Va, b,c £ M^. 

Letting / = <I> in the above expression yields 



{r xy) ■ {{ei x y) x {ly x v)) = J {r x {y x u)^)- . 

dn 



(2.52) 



an 



Gathering together (jOej), (jOSj), ([2^, p3T| )), and p32| ) yields 

3 



(i/ X y)(pi + ri 1 {v X y)(pi^ + / (^^ x y)'4'j 

af7 an ■?'=^ an 

m|3^ „ m„ I 

+ (^^ / ^"^^^ ^ 2/)Xj + n / (V(^i X y)xj) + ^Wp {Vipp X y)xj > 

J=i l*=i an an p=i an J 

1*=^ an an an J 

^ (li {yx + n {yx + ^'Wp {y x v)il)p \ 

an an an J 



r X Jqt. 



(2.53) 
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Combining (2.45) and (2.53) with the definitions of the matrices in (2.23)-( 2.30[ ), we obtain 



niQl 



Jqt- 



-Mi- Nr -C^'w - Wp{Lfl + R^^r + W^^w] 

l<p<m, 

-r X {Ml + Nr + C^'w) - rriQr x I, 

-N*i- Jr-C-'w- Wp{Lp + Rlr + W^w] 

l<p<m 

-I X {Ml + Nr + C^'^w) 

-r X {N*l + Jr + C'^w) - r x Jor. 



(2.54) 



This completes the proof of Proposition 2.3 



(2.55) 

□ 



2.4. Equations for the position and attitude. Now, we look at the dynamics of the position 
and attitude of the rigid body. We shall use unit quaternions. (We refer the reader to the 
Appendix for the notations and definitions used in what follows.) From (1.7) and (1.16), we 
obtain 

Q' = S{Qr)Q = QS{r), (2.56) 

with Q(0) = Id. 

Assuming that Q{t) is associated with a unit quaternion q{t), i.e. Q{t) = R{q{t)), then the 
dynamics of q are given by 

. 1 



-q * r 
2^ 



(see e.g. [S])- Expanding q as q = qo + q = qo + qii + q2j + 93^! this yields 



and 





qo + q 


■=2<- 


-q- r + qor 


+ 9 X 


r) 




( %] 




/ 90 


-91 -92 


-93 \ 




( ° ^ 




1 


91 


90 -93 


92 




ri 


92 


~ 2 


92 


93 90 


-91 






V ^3 J 




\ 93 


-92 91 


90 / 




\ ^3 / 



From (1.15), we see that the dynamics of h are given by 

h{t) = Q{t)l{t). 



Again, if Q{t) = R{q{t)), then (2.60) may be written as 

h = q * I * q* . 

Expanding q as q = qo + q = qo + qii + 52 J +93^, we obtain 

h= {qo + q)*l*{qo-q) =qo^ + 29o9 x / + (Z • q)q- qxl x q. 

and 



(2.57) 
(2.58) 

(2.59) 

(2.60) 
(2.61) 







h2 









7o + 9? 
2(9291 + 9093) 



92 - 93 2(gig2 - 9093) 2(^1^3 + 5092) 
9o - 9i + 92 - 93 2(^293 - 9091) 
2(9391 - 9092) 2(^3^2 + 9o9i) 9o - 9i - 92 + 93 




(2.62) 
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For q £ S"^, q may be parameterized by q, and it is thus sufficient to consider the dynamics 
of q which read 



IqW^r + q x r). 



The dynamics of h are then given by 



h={l- \\q\\'^)l + 2^/1 - ||g||Vx l + q)q- qxlxq. 
(Alternatively, one can substitute \/i — {qf + Q2 ^ ^i) ^'^ 9o hi both (2.59) and (2.62) 



(2.63) 
(2.64) 



2.5. Control system for the underwater vehicule. Using (2.57), (2.61), and Proposition 
we arrive to 

r h' = q*l*q*, 



2.3 



-q * r. 



(2.65) 



where {h, q, I, r, w) £ x x x R"^ x 



F{1, r, w) 



Sir) 
S{1) Sir) 



and 



C 



8-HCw' + Fil,r,w)), 



Cw] - Wp 

p=l 



(2.66) 
(2.67) 



For q £ S^ (i.e. Q £ 0), one can replace the two first equations in (2.65) by (2.64) and (2.63) 
respectively. This results in the system 



h' 



q 



(1 - II^IP)/ + 2v^l - II^IIVx / + (/ • q)q- qxlxq, 



r + g X r) 



(2.68) 



d-HCw' + Fil,r,w)). 



3. Control properties of the underwater vehicle 

3.1. Linearization at the equilibrium. When investigating the local controllability of a non- 
linear system around an equilibrium point, it is natural to look first at its linearization at the 
equilibrium point. 



To linearize the system (2.65 ) at the equilibrium point (/i, q, I, r, w) = (0, 1,0, 0, 0), we use the 



parameterization of Si by q, and consider instead the system (2.68). 
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h' 
2q' 



I, 



Proposition 3.1. The linearized system (3.1) with control w' € 
if, rank{C) = 6. 



(3.1) 



is controllable if, and only 



Proof. The proof fohows at once from Kalman rank condition, since {h,2q,l,r) G and 



rank 







Id 








2rank(C). 



□ 



Remark 3.2. It is easy to see that the controllability of the linearized system Iji3.1\) i mplies the 
(local) controllability of the full system (2.68). The main drawback of Proposition 3.1 is that the 
controllability of the linearized system (3.1) requires at least 6 control inputs (m>6). 

3.2. Simplications of the model resulting from symmetries. Now we are concerned with 
the local controllability of (2.68) with less than 6 controls inputs. To derive tractable geometric 
conditions, we consider rigid bodies with symmetries. Let us introduce the operators Si{y) = 
y - 2yiei for i = 1,2,3, i.e. 



S2{y) = (2/1,-2/2,2/3), 
Ssiy) = (2/1,2/2, -2/3)- 



(3.2) 



Definition 3.3. Let i £ {1, 2, 3}. We say that 0, is symmetric with respect to the plane {yi = 0} 
if 5i($7) = Q. Let / : J7 C — > M. // f{Si{y)) = e^jfiy) for any y £ 0. and some number 
ffj G { — 1, 1}; then f is said to be even (resp. odd) with respect to Si i/e^ = 1 (resp. = —1 ). 

The following proposition gather several useful properties of the symmetries Si, whose proofs 
are left to the reader. 6ip denotes the Kronecker symbol, i.e. 5jp = 1 if i = p, 6ip = otherwise. 

Proposition 3.4. Let i e {1,2,3}. Then 

(1) SiSi{a) = a,'iae M^" 

(2) Si{a) ■ Si{b) = a-b,ya,be R^; 

(3) Siia) X Si{b) = -Si{a x b), \/a,b £ R^. 

(4) IfSiin) = n, then i^iSiiy)) = Siii^iy)), Vy G dQ; 

(5) IffiSiiy)) = ef{y) with e € {±1}, then f{Si{y))v{Si{y)) = £Si{f{y)u{y)), My G 8^; 

(6) IfSi{n) = n, then Si{y) x iy{S,{y)) = -S,{y x u{y)), My G 8^1; 
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(7) Assume that S'j(il) = fl, and assume given a function g : dVl — t- M with g{Si{y)) = £g{y) 
for all y € dVt, where e £ {±1}- Then the solution f to the system 



in Q, 
g, on dQ, 



A/ = 0, 

dl 
dv 

^f{y) 0, as \y\ oo, 

which is defined up to an additive constant C , fulfills for a convenient choice of C 

f{S,{y)) = ef{y), ^y e n, 

V/(5,(y)) = eS,{Vfiy)), \/y € Q. 

(8) Let f and g be any functions that are even or odd with respect to Sp for some p G {1,2, 3}, 
and let h{y) = f{y)dug{y) ■ Then 



h{Sp{y))=ey^h{y) 



(3.3) 



p p p 



(9) Let f and g be as in (8), and let hi{y) = dif{y)dyg{y), where i G {1,2,3}. Then 

h,{Sp{y)) = {-l)''r'ePePhM, (3-4) 

(10) Let f and g be as in (8), and let hi{y) = (y x f{y))idyg{y), where i £ {1,2,3}. Then 

hi{Sp{y)) = -(-l)'^rePePh,iy), (3.5) 



I.e. £' 



V 



Applying Proposition 3.4 to the solutions (j)i, ipi, i = 1,2, 3, of (2.20)-(2.22), we obtain at once 
the following result. 

Corollary 3.5. Assume that Q is symmetric with respect to the plane {yp = 0} (i.e. Sp{^}) = Q) 
for some p G {1, 2, 3}. Then for any j G {1, 2, 3} 



(l^jiSpiy)) 



i.e. e^, = {—lypi, and 



V'jiSpiy)) 



4>j{y) if j T^p, 
-i)'"0i(y), 



-y^jiy) if 3 7^p, 

"Pjiy) if 3 =p, 
(-l)^-(^,(y). 



(3.6) 
(3.7) 

(3.8) 
(3.9) 



I.e. e 



-ly 



The following result shows how to exploit the symmetries of the rigid body and of the control 
inputs to simplify the matrices in ( 2.23 )-( 2.30) 
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Proposition 3.6. Assume that 0, is symmetric with respect to the plane {r/p = 0} for some 
p £ {1,2,3}. Then 



(1 
(2 
(3 
(4 

(5 

(6 

(7 



(9 
(10 
(11 



M,,=0 i/eiel = -1, I.e. 



Sip + Sjp = 1 (mod 2) ; 

ij — u y i^ipit^tpj ~ "I » " 



Jij — if Sip^s^j — 1, i.e. 



Nij = ifel el^ = -1, i.e. 



Sip + Sjp = 1 (mod 2); 

ij — >J y ~ " " 

Sip + Sjp = (mod 2); 

ij — ^ y ^(/,,^Xj ~ ^ 

(-i)^- = -4.; 

(C-^)ij = ife%^e% = -I, i.e. 

(-1^ = 4.; 



q Jij - ^ y l--^; ■"^'^'j^Xq ~ ^ " " 



(i^f),,■ = 0^/(-l)^-4,e5 =-1, ^.e. 



(T^f).,=0./(-l)^-4,.5, = -l, z. 



e. 



9-- - - - 

(i?,^),,- = 0z/ -(-1)^-6^,4, = -1, z.e. 



q 

1^<5^p ^ pP 



'Xj X? 

w/iere the matrices M, J, N, C*^, C-^, , , Wg*^ L^, and are defined in ([Z23])-| 
From now on, we assume that is invariant under the operators 5*2 and 5*3, i.e. 

Sp{n) = n, ype {2,3}, 



and that e^^ = 1, i.e. 



xiiSpiy)) = xiiy) yyedn,ypG{2,3}. 
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In other words, the set S and the control xi ^'^^ symmetric with respect to the two planes 
{y2 = 0} and {ys = 0}. As a consequence, several coefficients in the matrices in ( 2.23 )-( 2.30) 
vanish. 

More precisely, using (3.21 )-( 3.22 ) and Proposition 3.6, we see immediately that the matrices 
in (2.33) can be written 



M 




J 




(3.23) 



L 



M 



(^f)ll 








N 





{^1^)22 




A^32 








11 









33 









{Rr) 



32 




(3.24) 



(3.25) 



(3.26) 



and 



(H^f)ei 














[I 


{Li)23 








(Liu 



Ri 



(rOu 








{R{)22 








(^/)33 



(W^/)ei 




(3.27) 



(3.28) 



3.3. Toy problem. Before investigating the full system (2.68), it is very important to look at 
the simplest situation for which /ij = /j = for i = 2,3, q = 0, r = 0, and Wj = for j = 2, .., m. 



Lemma 3.7. Assume that (3.21 )-( 3.22 ) hold, and assume given some functions hi,li,wi G 
C"^([0,T]) satisfying 

' -- h 

(3.29) 

= aw'i+ filiwi+'yiwif, 



where 



a :- 



1 


\ K 




[ I'. 




)ii 



(^f)ll 



/3:= 



and 7 := 



-{Wi 



}ll 



mo + Mil mo + Mil mo + Mn 

Let h ■= {hi,0,0), q := (0,0,0), / := (/i,0,0), r := (0,0,0), and w := (i^i, 0, 0). Then 
{h, q,l,r,w) solves (2.68). 
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Proof. Let us set h = /iici, q = 0, I = hei, r = and w = {wi,0, ...,0), where {hi,li,wi) fulfills 
dOol ). From ( [3^ -( |335| ), we have that 

/ mo + Mil \ 








and 



This yields 



hdei = h 



Cw = wiCei 



-Wl 












Cw 



0. 



Replacing in (2.66), we obtain 



(3.30) 



(3.31) 



(3.32) 



F{l,r,w) 



p=i 



Lil + Rir + W:Iw 



-Wl 



-wi[ h 



Lfl + Wt'w 



Lfl + Wi'w 



J„ 



( (if )ii \ 








+ Wl 











(3.33) 



We conclude that (/i, g, /, r, w) is a solution of (2.68). 



□ 



Remark 3.8. If j + af3 = 0, then it follows from \\2\ Lemma 2.3] that for any T > we may 

associate with any pair {h\,h^) in a control input wi € C^(0,T) such that the solution 
{hi{t),li{t)) of (3.29) emanating from {h^,0) at t = reaches {hj,0) att = T. 
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3.4. Return method. The main result in this section (see below Theorem 3.11) is derived in 
following a strategy developed in [12] and inspired in part from Coron's return method. We first 



construct a (non trivial) loop-shaped trajectory of the control system (2.68), which is based on 
the computations performed in Lemma 3.7 (For this simple control system, we can require that 
wi{0) = 0, but we cannot in general require that wi{T) = 0.) Next, we compute the linearized 
system along the above reference trajectory. We use a controllability test from [T2| to investigate 
the controllability of the linearized system, in which the control appears with its time derivative. 
Finally, we derive the (local) controllability of the nonlinear system by a standard linearization 
argument. 

3.4.1. Construction of a loop-shaped trajectory. The construction differs slightly from those in 
|12j : indeed, to simplify the computations, we impose here that all the derivatives of Ii of order 
larger than two vanish at t = T. For given T > 0, let ^ G C°°(M; [0, 1]) be a function such that 

T 

if t < 

m 



3' 
2T 



ii t< 

1 if t > 

Pick any Aq > and let AG [— Aq, Aq] with A / 0. Set 

Mt) = m)it-Tf, iiit) = h[it), 

Note that 



t G 



/ii(0) = /ii(T)=Ii(0)=Ii(T) = 0, 
l[{T) = 2A / 0, if V) = for A; > 2. 
Next, define wi as the solution to the Cauchy problem 

^Z;i = a^^(li — fiiwi — ^wf), 
wi{0) = 0. 



(3.34) 

(3.35) 
(3.36) 



(3.37) 
(3.38) 



By a classical result on the continuous dependence of solutions of ODE's with respect to a 
parameter, we have that the solution wi of (3.37)-( |3^ ) is defined on [0,T] provided that Aq is 
small enough. Set ~h = (hi, 0,0), (0,0,0), w = (w i, 0, 0), / = (/i,0,0) and f = (0,0,0). 
According to Lemma 



3.7 



{h,q,l,r,w) is a solution of (2.68), which satisfies 
(h,ll,f){0) = = (h,llf){T). 

3.4.2. Linearization along the reference trajectory. Writing 

h 



q 

I 
r 



h + h, 
q + q, 

i + i, 

f + r. 



(3.39) 
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expanding in (2.68) in keeping only the first order terms in h,q, l and r, we obtain the following 
linear system 

h! = l + 2^x I, 
^1 



2' ' 



r ^ A{t) 



+ B{t)w + Cw'] , 



where the matrices A{t) G and B{t) G 



/ 

6x6 H^+^ ^ Tn,6xm defined as 



m 

B{t) 



dF 
~dl 

dF 

dw 



Setting 



we can rewrite (3.40) as 



h' 



P' 



m,rit)Mt)) 

m,f{t)Mt))- 

P = 2g, 

= I — I X p, 
= r, 

= d-'(A{t) 
ar 
k 



dF 
dr 



m,rit)Mt)) 



+ B{t)w + Cw' 



Obviously, (3.40) is controllable on [0, T] if, and only if, (3.44) is. Letting 

p J ' \ f 

we obtain the following control system 



/ 



We find that 



D{t) Id 
3-^A{t) 



+ 





3-'Bit) 



f + 







= : Ait) 



+ ^{t)f + ef. 



D 



-S(l) 




with S{1) 




B 





S{1) 



C-wi 



J M „ T M „ 

- wi 



f 



(3.40) 



(3.41) 
(3.42) 

(3.43) 
(3.44) 



(3.45) 







-h , 

/ 








T J P-, 1 






Wiei 







A 



(^f)llwl 







































As 






















As 





iVs2ll - {Ri)22Wl 





A&2 
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and that 

A2Q 






V A2 -iVssIl - (^f )S3W^1 / 

with 

Ae = -(mo + Mn)IT- ((Oii + (iZf)2s)^I'i, 
A5 = (mo + Afn)IT+((C*0ii-(<)s2)u'i, 
As = (M33-Mn)II- ((C*^)n + (Lf)23)Wi, 

A2 = (Mn-M22)IT+((C^')ll-(Lf)32)Wi. 

From now on, we suppose in addition to ( 3.21[ )-( p.22[ ) that xi is chosen so that 

a / 0. (3.46) 

3.4.3. Linear control systems with one derivative in the control. Let us consider any hnear con- 
trol system of the form 

x = A{t)x + 'B{t)u + Cu (3.47) 

where E M'' is the state (n > 1), n e is the control input (m > 1), ^1 G C°^([0, T]; M"^"), 
S E C°°([0,r];]R"^'"), and 6 G M"''™. Define a sequence of matrices Mi(t) G M"^"" by 

Mo(0 = B(t) +yi(t)e, and JAi{t) = Mi_i(t) - yi(t)Mi_i(t), Vi > 1, Vt G [0,T]. (3.48) 

Introduce the reachable set 

3^«(o)=o = {a^T G M"; 3u G i?^(0,r;M™) with u(0) = such that 



XT = 3^(7"), where x(-) solves (3.47) and x(0) = 0}. 

Then the following result holds. 

Proposition 3.9. |Ll2i Propositions 2.4 and 2.5] Let e > A e C"^((- £,r + e);R"^") and 
"B G C"^((-e,r + e);M"^"^), and let {Mi)i>o be the sequence defined in ( |3.48D . Then for all 
to G [0,r], we have that 

= eM" + S'pan{0(r, to) Mi{to)u; n G M™, i > 0}, (3.49) 
where (p denotes the fundamental solution associated with the system x = A{t)x. 

Recall that the fundamental solution associated with x = A{t)x is defined as the solution to 

^ = A{t)cl){t,s), 
ch(s,s) = Id. 



CONTROL OF UNDERWATER VEHICLES IN INVISCID FLUIDS. 



I: IRROTATIONAL FLOWS. 23 



For notational convenience, we introduce the matrices 

A{t) = 3'^A{t), B{t) = 3-^B{t), C = 3-^C, 



Viit) 



where A{t) £ 



while 



»6x6 



Bit),C,Uiit),Vi{t) G 

Uo{t) 
Vo{t) 



»6xm 



Then 



C 

I3{t) + A{t)C 



(3.50) 



(3.51) 



Ui{t) UUit) - Dim-lit) - Vi-iit) 

v^it) J \ y/„i(t)-i(t)y,_i(t) 

In certain situations, half of the terms Uiit) and Viit) vanish at t = 
whose proof is given in Appendix, will be used thereafter. 



(3.52) 

T. The following result, 



Proposition 3.10. If C £ M^^™' is given and A, D (resp. B) denote some functions in 
C°^([0,T];]R6x6) ^^g<,p_ in C^^ilO, T]-R^'''^) fulfilling 



then the sequences (f7j)i>o and (Fi)i>o defined in (3.51 )-(3.52) satisfy 



vif\T)- 



v^I'^'\t) 



'2k+l 
'2k 



0, 



0, 



VA;,/ G N, 
VA;,/ G N. 



(3.53) 

(3.54) 
(3.55) 



The following result, which is one of the main results in this paper, shows that under suitable 
assumptions the local controllability of (2.68) holds with less than six control inputs. 



Theorem 3.11. Assume that ( |3.2l[ ), ( |3.22[ ) and ( |3.46[ ) hold. Pick any T > 0. // the rank 
condition 

rank (C, Mo(r), Mi(r), M2(r), ...) = 12 (3.56) 

holds, then the system (2.68) with state (/i, q, Z, r) G M"^^ and control w G M™ is locally controllable 
around the origin in time T . We can also impose that the control input w G -ff^(0, T; M™) satisfies 
If (0) = 0. Moreover, for some r] > 0, there is a C map from S]g24(0,7?) to i/2(o,r;M™), which 
associates with ih^, qo,lQ,rQ, hr, qTjT,rT) a control satisfying wiO) = and steering the state 
of the .system from (/iq, qo, ^o, tq) at t = to (Ht, qr, It, ^t) at t = T . 

Proof. Step 1: Controllability of the linearized system. 

Letting to = ^ in Proposition 3.9 yields 



/(o)=o 



+ Y,UiiT)l 

i>0 



t>12 



Thus, if the condition (3.56) is fulfilled, we infer that 3?j(o)=o 
controllable. The same is true for (3.40). 
Step 2: Local controllability of the nonlinear system. 

Let us introduce the Hilbert space 

i^^ X {/ G H^iO,T;R"'); /(O) = 0} 



i.e. the system (3.44) is 
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endowed with its natural Hilbertian norm 

Wix, f)\\li = \\x\\li2 + ||/|Ih2(o,t) • 

We denote by i3j{(0, 6) the open ball in with center and radius 6, i.e. 

Bji{0,6) = {ix,f)eJi;\\ix,f)\\:^<6}. 

Let us introduce the map 

i{ho,qo,lo,ro),f) ^ {ho,qo,lo,ro,h{T),q(T),l(T),r{T)), 
where {h{t),q{t),l{t),r{t)) denotes the solution of 

h' = (1 - ||g||2)/ + 2Y^l - \\q\\^ q x I + {I ■ q)q-qx I X q, 

q = - ||g IP r + g X r), 



(3.57) 

r\C{w' + f) + F{l,T,w + f)), 



[ (/i(0),g(0),/(0),r(0)) = (/io,«),/o,r-o). 

Note that T is well defined for 5 > small enough (provided that Aq has been taken sufficiently 
smah). Using the Sobolev embedding H'^{Q,T;W) C C1([0, T]; M™), we can prove as in [231 
Theorem 1] that T is of class on i?j{(0,(5) and that its tangent linear map at the origin is 
given by 

dr(0)((/io,?o,/o,ro),/) = {ho,?o,io,ro,HT),?iT)JiT),fiT)), 
where {h{t), q{t), l{t),r{t)) solves the system (3.40) with the initial conditions 

(/i(0),?(0),/(0),f(0)) = (/io,?o,Zo,ro). 



We know from Step 2 that (3.40) is controllable, so that cir(O) is onto. Let V := (ker cir(O)) 
denote the orthogonal complement of ker dT(0) in IK. Then (ir(0)|y is invertible, and therefore it 
follows from the inverse function theorem that the map T\v '■ V ^ M?^ is locally invertible at the 
origin. More precisely, there exists a number S > and an open set cj C M^^ containing 0, such 
that the map T : -Bj{(0, 6)r\V — )• w is well-defined, of class C^, invertible, and with an inverse map 
of class C^. Let us denote this inverse map by F"^, and let us write r~^(xo, xt) = {xq, /(xq, xt))- 
Finally, let us set w = w + f . Then, for > small enough, we have that 

w G C7i(SiR24(0,r/),F2(0,T;M'")), (3.58) 

and that for || (/io; ^O; ^O) '''Oi ^Ti ^Ti ^Tj '"t)! |r24 < ry, the solution {h{t),q{t),l{t),r{t)) of system 
(2.68), with the initial conditions 

(/i(0),g-'(0),/(0),r(0)) = iho,qo,lo,ro), 

satisfies 

{h{T),q{T),l{T),r{T)) = {hT,qT,lT,rT). 



The proof of Theorem 3.11| is complete. □ 
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We now derive two corollaries of Theorem 3.11, that will be used in the next section. We 
introduce the matrices 



Al 



Ar 



where 



Al 



( -{^\ 











a(M-_ 



11 



12) 









-((^/)32-(C^%: 



(3.59) 






-(i^f)33 



a(M33 - Mn) - {{L{)2z + (C^^ii) 




Ar 











11 

aiV32 - {R{)22 











B 



-a5(ei)C 



M 



and 



W(ei 



D 



ei 
Liei 



Liei 



\ 

-«)23 






-aiV23 - (^/)33 ) 



(3.60) 



(3.61) 



M 



(3.62) 








(3.63) 



The first corollary will be used later to derive a controllability result with only four control 
inputs. 



(3.64) 



(3.65) 



Corollary 3.12. // both rank conditions 

rank {C,B + Ad'^C) = 6 

and 

rank (C, ^d'Dd'^C + B + A^^^C) = 6 

are fulfilled, then the condition ( 3.56| ) is satisfied for any T > 0, so that the conclusion of 
Theorem 3.11 is valid for any T > 0. 

Proof. We distinguish two cases. 
Case 1: 7 + a/3 = 0. 

We begin with the "simplest" case when 7 + q/3 = 0. Pick any T > and let Ii, be as in 
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(3.34) and ( 3.37 )-( 3.38 ). Let := li — awi. It is clear that = (3wigi, hence gi = 0. We infer 
that 

w[^\t) = a-n[^\T) = for A:GN\{1}, 
w[{T) = a-%{T) = 2A/a / 0. 

It follows that 

^W(r) = 0, S(^-)(r) = 0, D('=)(r) = forA;GN\{l}, (3.66) 

A!{T) = w'^{T)A, B'{T) = w[{T)B, D' (T) = w[{T)U . (3.67) 

Applying Proposition |3.10] we infer that 
rank (S, Mo(r), Mi(r), M2(r)) 

= rank ( ( ^ ) , ( ^ ^^r^^^j ) , ( v,It) 

On the other hand, it is easily seen that 

Vi{T) = V^{T) = r^B'{T) + r^A'{T)r^C = w[{T) + 3'^A3-^C) , 

U2{T) = -D'{T)Uq{T) - 2V^{T) = -w\{T)[V^d'^C + 20-1 (B + Ag-^C)]. 

It follows that 

rank(C', Vi{T)) = rank (C, B + Ad'^C) = 6, 
rank(C', U2{T)) = rank (C, ^dm^^C + B + AJ'^C) = 6, 



U2{T) 




and 



rank (C, Mo(r), Mi(r), M2(r)) = 12. 



Thus (3.56) is satisfied, as desired. 

Case 2.7 + a/3 7^ 0. We claim that for T > arbitrary chosen and Aq small enough, we have 
for < A < Ao, 

rank (C, Mo(r), Mi(r), M2(T)) = 12. 



First, \\h\\w^,^io,T) = 0{X) still with li{T) = li{T) = 0. From ( |3.37D -( |3.38D , we infer with 
Gronwall lemma (for Aq small enough) that wi is well defined on [0,T] and that ||^^i ||l°°(o ,t) = 
0(A). This also yields (with ( |3.37| )) | |tt;i| |v(/2,oo(o = 0(A). Next, integrating in (3.37) over 
(0,r) yields wi{T) = 0{X^). Finally, derivating in ( |3.37D gives %i{T) = 0{X^). We conclude 
that 



while 



(A(r), B{T),A{T), B{T)) = 0(A2), D{T) = 0, 
(i(r), i?(r), DiT)) = (2A/a)(A, B, D) + O(A^), 
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for Ii(r) = aWi(T) + 0(A4). It follows that 
rank (C, Mo(T), Mi(r), M2(T)) 



rank 













r M3Dr + 2(B + Ar ^c)] 





12, 



for < A < Ao with Aq small enough, as desired. □ 

The second one is based on the explicit computations of Mj(T) for i < 8. It will be used later 
to derive a controllability result with only three controls inputs. 

Corollary 3.13. Let E := B + A^'^C. // the conditions 

rank (C, E, Aa^^E, {Ar^f'^', (Aa~i)3E) = 6, (3.68) 

and 

rank (C, ]^dm'^C + E, {md'^ + 2Aa-i)E, 

(saoa^^ + iiAa"^)Aa"^E, (iraoa^^ + 64Aa~^)(Aa"^)^E) = e, (s.eo) 



are fulfilled, then the condition (3.56) is satisfied, so that the conclusion of Theorem 3.11 
valid. 



IS 



Proof. The proof is almost the same as those of Corollary 3.12, the only difference being that 
we need now to compute Mj(T) for i < 8. In view of Proposition 3.10, it is sufficient in Case 1 
(7 + a/3 = 0) to compute ^^(r) for i £ {1, 3, 5, 7} and Ui{T) for i £ {2, 4, 6, 8}. The results are 
displayed in two propositions, whose proofs are given in Appendix. 



Proposition 3.14. Assume that the pair {hi,li) is as in (3.34), that wi is as in (3.37)-(3.38), 
and that 7 + a/3 = 0. Then we have 

-M'{T)V^{T), 



ViiT) = V^{T), (3.70) 

V3(,T) = -3A'{T)V^{T), (3.71) 

V^iT) = 15i'(T)Vo'(r), (3.72) 

Vr{T) = -105A'{TfVQ{T). (3.73) 

Proposition 3.15. Assume that the pair (hiji) is as in (3.34), thatwi is as in (3.37) -( [338 ), 

and that that 7 + a/3 = 0. Then we have 

U2{T) = -D'{T)Uo{T)-2V^{T), (3.74) 

Ui(T) = A{D'{T) + 2A'{T))Vq(T), (3.75) 

Uq{T) = -3{8D'{T) + 11A'{T))A'{T)Vq{T), (3.76) 

UsiT) = 6(i7L>'(r) + 64i'(r))i'(r)Vo(r). (3.77) 
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□ 



4. Examples 

This section is devoted to examples of vehicles with "quite simple" shapes, for which the 



coefficients in the matrices in ( 2.23 )-( 2.30) can be computed explicitly. We begin with the case 



of a vehicle with one axis of revolution, for which the controllability fails for any choice of the 
flow controls. 

4.1. Solid of revolution. Let / € ei([a,6];M) be a nonnegative function such that /(a) = 
f{b) = 0, and let 

S= cos sine); yi € [a,b], sG [0,1], [0,27r]}. 

In other words, S is a solid of revolution (see Figure [T]). 

Assume that the density p depends on yi only, i.e. p = p{yi). Clearly Jq = diag(Ji, J2, J2)- 
On the other hand, 

5f]= {(yi,/(yi)cos0,/(yi)sine); yi G [a, 6], ^G [0,27r]}, 

and the normal vector v to dVt is given by 

1 



so that 



i'{yi,G) = ^==^== /'(yi),-cos6', -sin6' 

Vl + (/ (yi)) ^ 

{y X u){yi,e) = — ===^fO,sm0,-cos0 



It follows that {y X v) ■ ei = 0. Replacing in (1.24), we obtain 



Jiri = [Jq-t) ■ ei = -(r x Jor) • d = J2r2r3 - J2r2rz = 0, 
which indicates that the angular velocity ri is not controllable. 




Figure 1. A solid of revolution. 
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4.2. Ellipsoidal vehicle. We assume here that the vehicle fills the ellipsoid 

§={ye M^; (yi/ci)2 + (y^/c^f + (^3/03)2 < 1} (4.1) 

where ci > C2 > C3 > denote some numbers. Our first aim is to compute explicitly the 



functions (pi and ifi for i = 1,2,3, which solve (2.20)- p.22[ ) for 

n={ye M3; (yi/ci)2 + (^2/02)^ + (2/3/03)' > !}• 



4.2.1. Computations of the functions (pi and (fii. We follow closely [16l pp. 148-155]. We introduce 
a special system of orthogonal curvilinear coordinates, denoted by (A, /i, i^), which are defined 
as the roots of the equation 



4 + ' 



+ 



+ 



vl 



1 = 



viewed as a cubic in Q. It is clear that (4.2) has three real roots: A € (— C3, +00), /U G 
and V G (— cf , — Cj). 



It follows immediately from the above definition of A, /i, i^, that 



vl -I 

Cl + e (c2 + 0)(c2 + e)(ci + , 



This yields 



vl 



vl 



vl 



{cl + \){cl + ^l){c\ + v) 

{cl-cl){cl-cl) ■ 

{cl + \){cl + ^JL){cl + v) 

{c\-cl){cl- cl) ■ 

{cl + \){cl + lM){cl + v) 



{cl-cl){cl-cl) 



d\yi 
d\y2 
dxys 



yi 



2{ci + X) 
1 y2 

2{4 + xy 

1 ys 

2{ci + xy 



We introduce the scale factors 



(A-^)(A-z.) 



(A + c2)(A + ci)(A + ci)' 



(4.2) 



-2) 



-3^ 



(4.3) 



1 / (^-„){^-X) 

2^ {^, + cl){^, + 4){^, + 4y 



1 / („-X){u-fj.) 



and the function 



/(A) = ^(A + c2)(A + ci)(A + ci). 
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If ^ is any smooth function of A, then its Laplacian is given by 

according to [161 (7) p- 150]. We search cpi in the form (/>i(?/i, 2/2, ys) = 2/2, ys)- Then 

= A4>i = yiA(i + 2di^i. 
Assuming furthermore that depends only on A, we obtain that 



Vi h\ cj + Xh\ c2 + A (A-^)(A-i/)' 



Combining (4.6) with (4.5) and (4.7), we arrive to 



= dx{f{\)dxii) + -^f{\)^X^^, 



which is readily integrated as 



+ 00 



-a 



ds 

{cj + s)f{s) 



+ C. 



(4.5) 
1 

(4.6) 
(4.7) 



A = 0, then K2\\ reads 



We choose the constant C = for (2.22) to be fulfilled. As dil is represented by the equation 

OuCPi = ti ^ Oxt,i = . 

Vi Vi 



We infer that Cj = ciC2C3/(2 — Oj), where 



tti = C1C2C3 



ds 



It is easy seen that 



2C2C3 2ciC2 

2^ < Cti < 

3c-i 3cq 



-1 

It follows that if ci,C2,C3 are sufficiently close, then Oi is different from 2, so that d is well 
defined. We conclude that 



2 - a, 



-yi, Vy € (917. 



(4.8) 



Let us now proceed to the computation of ipi. We search ipi in the form ipi{y) = ^^^^^ ^^(y), 
where depends only on A. We obtain 

A^, + 2 ^ ^ = 04^9A(/(A)aA^.) + ( 5] ^^)/(A)aA^, = 0, 



and hence 



+00 



ds. 
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Prom ([2^- ([2^, we infer that 

Cl 



Co 



C. 



C1C2C3 


[4 - 4) 


2 - 


-/3i 


C1C2C3 


[4 - 4) 


2 - 


-h 


C1C2C3 


[4 - 4) 


2 - 





+00 



^1 = C1C2C3(C2 + C3) 



f32 = C1C2C3(C3 + Cl) 



ds 





+00 



{4 + s){ci + s)f{s) 





+00 



133 = CiC2C3{cl + C2) 







+ + ,)/(,) ^ 

ds 

{cl + s){cl + s)f{sy 



Note that at the limit case ci = C2 = 03, we obtain Pi = (32 = (3^ = 4/5. Therefore, if ci,C2 
and C3 are near but different, then /3j is different from 2, and therefore Ci is weh defined. We 
conclude that 

+00 

yiy2y3 







cf + s 
Pis) 



ds 



Vi 



(4.9) 



4.2.2. Controllability of the ellipsoid with six controls. Assume still that S is given by (4.1 ). Note 
that S is symmetric with respect to the plane {up = 0} for p = 1,2,3. Assume given six functions 
Xj, j = 1, . . . , 6, each being symmetric with respect to the plane {up = 0} for p = 1, 2, 3, with 



'XI 



'X4 



-X2 



-X5 



'X3 



'X6 



1 


p = 


1 


1 


p = 


2 


-1 


p = 


3 


-1 


p = 


1 


-1 


p = 


2 


1 


p = 


3 



(4.10) 



To obtain this kind of controls in practice, we can proceed as follows: 

• We build six tunnels in the rigid body, as drawn in Figure [2} 

• We divide the six tunnels in three groups of two parallel tunnels; that is, we put together 
the tunnels located in the same plane (see Figure [s]). 

• Let wi and W2 denote the effective flow controls in the two tunnels located in the plane 
{?/3 = 0}. They may appear together in ( |1.21[ ) as ^2, 2/3) + mxivi, -y2,yz), 
where x ^ C°°{dQ) is some function with 

Supp X C {2/2 > 0}, = —1, and 

We introduce the (new) support functions 

xi(yi,y2,?/3) = x{yi,y2,yi) + x{yi,-y2,yz) 
X6(yi,2/2,y3) = x{yi,y2,y3) - x{yi,-y2,y3) 



1. 
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Figure 2. Ellipsoid with six controls. 




















■■. Ill 











Figure 3. Independent controls in each plane, 
and the (new) control inputs 



Wq 



Wl + W2 

2 ' 

Wl — W2 



(See Figure |4j) Then (4.10) is satisfied for xi and xe, and 

wixiyi,y2,y3) + W2xiyi,-y2,y3) = wixiiyi,y2,y3) + weX6iyi,y2,y3) 

The same can be done in the other planes {yi = 0} and {y2 = 0}. 
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Wi 



Wi + W2 



W2 



W2 



Wl - 'W2 



Figure 4. Definition of the new controls in the plane {y^ = 0}. 



We notice that C is a diagonal matrix: 

C = -diag(Ci,C2,C3,C4,C5,C6), 

with 

Ci= I (j)iXi, « = 1,2,3, and Cj+3 = J (piXi+3, « = 1,2,3. 



an 



Prom (|4.8|)-(|4.9|), there are some constants Cj 7^ 0, i = 1, ...,6, which depend only on ci,C2 and 



C3, such that 
Ci = Ci 



ViXiiy), C'i+3 = C'i+3 



(9f7n(0,+oo)3 



ar2n{o,+oo)3 



Xi+3(y), i = 1,2,3. (4.11) 



By (4.11), we have that Cj 7^ for i = 1, -..,6, and hence rank(C) = 6 if, in addition to (4.10) 
it holds 



Xi^O, « = 1, •••,6, 

Xi>Oondnn (0, +oo)^ i = 1, 6. 



(4.12) 
(4.13) 



By Proposition 3.1 and Theorem 13.11 it follows that both the linearized system (3.1) and the 



nonlinear system (2.68) are (locally) controllable. 



Remark 4.1. Since e^^ = (-1)''" = -1, we have that {Lf)ii = {W^)ii = 0, and hence 
/3 = 7 = 0. Thus J + a(3 = 0. Proceeding as in [12', Theorem 2.2], one can prove that, 
under certain rank conditions, two arbitrary states of the form (h, q, 0, 0) can be connected by 
trajectories of the ellipsoid in (sufficiently) large time. 
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In the following sections, we shall be concerned with the controllability of the ellipsoid with 
less controls (namely, 4 controls and 3 controls). If, in the pair (xiiXe)) only xe is available, 
then X6 can be generated as above by two propellers controlled in the same way (Figure [5] left), 
or by only one propeller by choosing an appropriate scheme for the tunnels (Figure [5] middle) . 
In what follows, to indicate that the flows in the two tunnels are linked, we draw a transversal 
line in bold between the two tunnels (Figure Is] right). 




Figure 5. Two ways to generate xe- 



4.2.3. Controllability of the ellipsoid with four controls. We consider the same controllers Xi^Xi-, Xb 
and X6 as above, still satisfying ( 4.10[ ), (4.12), (4.13). (See Figure [g]) If the density p is scaled 
by a factor A, i.e. p{x) is replaced by p^{x) = Xp{x) where A > 0, then the mass and the inertia 
matrix are scaled in the same way; that is, niQ and Jq are replaced by 

tA 



Xmn 



Jr. — AJn 



Control xi 



Control X4 







V-3 






























Control X5 



Control X6 

!/3 


















j5> 










-■•.!/2 





Figure 6. Ellipsoid with four controls. 
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Thus, if A — )■ oo, then uiq — t- oo, [Jq]'^ 0, and [3^]-'^ 0. (Note that large values of A 
are not compatible with the neutral buoyancy, but they prove to be useful to identify geometric 
configurations leading to controllability results with less than six control inputs.) 

Note that the matrices M, J, N, C^^ C-^ , L^^ W^'^ , L^,Rp, keep constant when A 
oo. In particular. 



lim Ad'^C = 0, 



1 



lim-aDa^^c 



Let S~ 



hm B. Then B°° 

A— >oo 



and C are given by 



0. 



/ 













\ 



\ 

Be 




J 



c 



( c, 












d 





V 







\ 











Ce J 



with 



55 



an 



l-VV'5)z^3, Be 



an 



Thus, if i^s 7^ and Be / 0, we see that (3.64) and (3.65) are fulfilled, so that the local 



controllability of (2.68) is ensured by Corollary 3.12 for A large enough. We note then that the 
matrix in M^^^ obtained by gathering together the four columns of C and the last two columns 
of -B°° is invertible. Let i?i(A) G M^^^ (resp. i?2(A) G M^^^) denote the matrix obtained by 
gathering together the four columns of C with the last two columns of B + A3~^C (resp. with 
the last two columns of \3T>3'^C + B + Ad'^C). Then for A > 1, we have 

deti?i(A)/0 and det i?2(A) / 0. 

Since the coefficients of i?i(A), i?2(A) are rational functions of A, we infer that the equation 

det Ri{\) -det i?2(A) = 

is an algebraic equation in A. Therefore, it has at most a finite set of roots in (0, +oo), that we 
denote by ^critical- We conclude that for any A G (0, +oo) \ AcrUicai, the local controllability of 



(2.68) still holds. In particular, we can consider values of A arbitrary close to the value A = 1 
imposed by (1.10). The issue whether 1 G ^-critical seems hard to address without computing 



numerically all the coefficients in our system. 



4.2.4. Controllability of the ellipsoid with three controls. Assume that Xi5X4)X5 ^tnd xe are as 
above (satisfying (4.10), (4.12), (4.13)), and consider now the controls supported by xi) X4 and 
X5 = X5 + X6 (see Figure |7|). 
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Control xi Control Xi 





Control X5 



Figure 7. Ellipsoid with three controls. 



Doing the same scaling for the density, and letting A 
C read 



oo, we see that the matrices and 



5= 











\ 








^6 










^5 


































/ 



c 



( c, 






V 



C4 




\ 





C5 
Ce ) 



where the coefficients B^, Bq,Ci,C4,C5,Cq are as above. For simplicity, we assume that the 
principal axes of inertia of the vehicule coincide with the axes of the ellipsoid. Then the matrix 
Jo is diagonal (see |3| ) with entries Ji, J2, Ja- Notice that the first and fourth coordinates are 
well controlled (using xi and X4)i and that the other coordinates are decoupled from them, at 
least asymptotically (i.e. when A — )• 00). Let = lim;^_i.oo A (i.e. is obtained by letting 
a = in A). Let K £ M^^*^ denote the matrix obtained from by removing the ffist and 
fourth lines (resp. columns), and let 6 G (resp. c G M'^) denote the vector obtained from the 
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last column of (resp. C) by removing the first and fourth coordinates, namely 



-{L 1^)22 



V -((L/)32-(C^^ 



K 







33 



/^6\ 



V / 













-iRi)22 





32 



-«)23 \ 




-(^/)33 J 



Let finally 



F 



1 ^ 
















1 

























/ 



and G 



/ 













\ 

























J2' 







V 











^3-^ 


/ 



Then, keeping only the leading terms as A — )• 00, we see that (3.68) holds if 

rank (c, 6, KGh, {KOfb, (KGfb) = 4 



(4.14) 



while (3.69) holds if 



rank ( c,b, [iC^^)uF + 2K]Gb, [8{C''")uF + llK]GKGb 



[17(C7^^)iiF + Q4:K]G{KGYb ] = 4. (4.15) 



Note that (4.14) is satisfied whenever 

rank(5, KGb, {KGfb, (KGfb) = 4, 



(4.16) 



which is nothing but the Kalman rank condition for the system x = KGx + bu. However, it 



is clear that we should take advantage of the presence c in (4.14). As previously, this gives a 
controllability result for A ^> 1, but such a result is also valid for all the positive A's except those 
in a finite set defined by an algebraic equation. 



5. Appendix 

5.1. Quaternions and rotations. Quaternions are a convenient tool for representing rotations 
of objects in three dimensions. For that reason, they are widely used in robotic, navigation, 
flight dynamics, etc. (See e.g. [HEl!)- We limit ourselves to introducing the few definitions and 
properties needed to deal with the dynamics of h and Q. (We refer the reader to [1] for more 
details.) 

The set of quaternions, denoted by H, is a noncommutative field containing C and which is a 
R-algebra of dimension 4. Any quaternion g G IHI may be written as 



q = qQ + qii + q2j + qsk, 
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where {qo, qi,q2-,Qz) G and i, j, A; G H are some quaternions whose products will be given later. 
We say that go (resp. qii + g2j + is the real part (resp. the imaginary part) of q. Identifying 
the imaginary part qii + g2j + Qsk with the vector q = ((?i,(?2,93) G we can represent the 
quaternion q as q = [qo, q], where go £ 1^ (resp. q £ M^) is the scalar part (resp. the vector part) 
of q. The addition, scalar multiplication and quaternion multiplication are defined respectively 

by 

\po,p] + [qo,q] = \po + qo,p + q], 
t[qo,q] = [tqo,tq], 

[po,p] * [qo,q] = [poqo-p- q,poq + qop + p x q], 

where "•" is the dot product and "x" is the cross product. We stress that the quaternion 
multiplication * is not commutative. Actually, we have that 

i* j = k, j * i = —k, 
j * k = i, k* j = —i, 
k* i = j, i* k = —j, 

Any pure scalar go and any pure vector g may be viewed as quaternions 

qo = [qo,o], q = [o,q], 

and hence any quaternion g = [go, g] can be written as the sum of a scalar and a vector, namely 

q = qo + q- 

The cross product of vectors extends to quaternions by setting 

px q = ^{p* q- q*p) = [o,px g]. 

The conjugate of a quaternion g = [go, g] is g* = [go, — The norm of g is 

Ikll = (kol^ + 

From 



we infer that 



q*q =q *q = \\q\\ , 



-1 q 

1 = TTZU2- 



\\q\ 

A unit quaternion is a quaternion of norm 1. The set of unit quaternions may be identified with 

S^. It is a group for *. 

Any unit quaternion g = [go, g] can be written in the form 

oi . a ^ . 
g = cos — + sm — u, (5.1) 

where a G M and u G with ll^uH = 1. Note that the writing is not unique: if the pair (a, « ) is 
convenient, the same is true for the pairs (—a, —H) and {a + Ak7r,il) (A; G Z), as well. However, 
if we impose that a G [0, 27r], then a is unique, and u is unique for |go| < 1. (However, any 
u e is convenient for |go| = 1.) 
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For any unit quaternion q, let the matrix R{q) £ M'^^^ be defined by 

R{q)v = q*v*q* Vv G (5.2) 

Then R{q) is found to be 

Qo + Qi-qI- I3 HqiQ2 - qoQs) 2(^1^3 + qoq2) 

= \ 2(9291 + ^0^3) qi-qj + qi-qi 2(^2^3-9091) 

2(9391-9092) 2(93^2+9091) 9o- 91-92+ 93 



For q given by (5.1), then R{q) is the rotation around the axis Mu of angle a. 
Note that R{qi * 92) = R{qi)R{q2) (i-e. i? is a group homomorphism) , hence 

R{l) = Id, R{q*) = R{q)-^ . 

We notice that the map 9 — t- R{q) from the unit quaternions set to 5*0(3) is onto, but not 
one-to-one, for R{—q) = R{q). It becomes one-to-one when restricted to the open set 

5"+ := {9 = [90,9] e H; Ikll = 1 and 90 > 0}. 

Furthermore, the map is a smooth invertible map from 5^^ onto an open neighbourhood of 
Id in SO{3). On the other hand, the map 

q^q = [90, 9 ] = - II9IP, 9 ] 

is a smooth invertible map from the unit ball -Bi(O) = {9 G M^; H^H < 1} onto S'^. Thus the 
rotations in can be parameterized by 9 G i?i(0). 

5.2. Proof of Proposition 3.10| , Let us prove by induction on /c G N that 

V^^^\t) = VZgN. (5.3) 

The property is clearly true for k = 0, since 

fJ^'^(T) = 5(2') (T) + A^'^^\t)c = 0, 

by (3.53). Assume that (5.3) is established for some /c G N. Then by (3.52) applied twice, we 
have 

V2k+2 = - 2il^2fc - ^'^2fc + iV2fc, 

llGnCG 

v!Sl^(T) = - 2(iy4)(20(r) - {a'V2,)('^Ht) + (iV2,)(20(r). (5.4) 

The first term in the r.h.s. of ( |5.4[ ) is null by ( |5.3| ). The second one is also null, for by Leibniz' 
rule 

21 

and i(P)(r) = if p is even, while ~^+^^(r) = if p is odd. One proves in a similar way 
that the third and fourth terms in the r.h.s. of (5.4) are null, noticing that for p odd we have 

(i2)(p)(T) = 2{AA'Yp-'^\T) = 0. (5.5) 



From (5.3), we infer that 
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Let US proceed to the proof of (3.55). Again, we first prove by induction on /c G N that 

V/ G N. 



It foUows from ( [3l5l| ), ( [3l52| ) and that 



(5.6) 



yJ2')(r) = 



V/t e N. 



Assume that (5.6) is true for some G N. Then, by (3.52) apphed twice, 
uSUt) = US!^P{T) - (DC/2fe+i)(2'+i)(r) - vi'^,'\T) - {DU2t+2)^'HT) - vi%(T). (5.7) 



Using ( 3.53|) , (3.54) and (5.6), we see that all the terms in the r.h.s. of ( |5.7[ ), except possibly 
(DC/2fe+^(T), are null. Finally, 



Using Leibniz' rule for each term, noticing that in each pair (p, q) with p + q = 21, p and q are 
simultaneously even or odd, and using (3.53), (3.54), (5.5) (with A replaced by D), and (5.6), 
we conclude that {DU2k+2)^'^^\T) = 0, so that U^i%iT) = 0. 

is obvious for k = 0, while for k > 1 

j{2l + 2)frj,^ Ij^jj^ .\{'^1 + '^){T) -V^I^^^\T) =Q 

is 



Finally, uf^^^\T) 



by (3.53), (3.54) and ( |5.6| ) (with 2k + 1 replaced by 2k — 1). The proof of Proposition 3.10 
complete. 

5.3. Proof of Proposition |3.14| , From ( |3.50| ), ( |3.51| ) and ( |3.66D , we obtain successively 
ViiT) = Vf(,{T) = B'{T) + A'{T)C = W{T){3~^B + 3~^A3~^C) 

{v; - ivi)'(r) 

{V^ - AVoTiT) - {AV^YiT) 
V^'{T) - 2A'{T)V^{T) - A'{T)Vi{T) 
-3A'{T)Vo{T). 



Successive applications of (3.51) yield 



Vj{T) 



1=0 
5 



(5. 



(5.9) 



i=0 



Since vi''\T) 



for k > 2, it remains to estimate the terms {AVi)'^'^-''\T) and {AViY^~'\T). 



Notice first that by (3.66) and Leibniz' rule 



{AViY''\T) = kA'{T)v}''~^\T). 
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Thus, from (3.66) and (3.71), we have that 



(5.10) 

(iyi)(3)(r) = 3A'{T)Vl'{T) = 3A'{T){vi^\T) - {AVo)"{T)) = -GA' {TfV^{T), (5.11) 



(^yo)^'^^) = 0, 



(5.12) 
(5.13) 

(5.14) 
(5.15) 
(5.16) 



{AV2)"iT) = 2A'{T)Vi{T) = 2A'{T)V3iT) = -6A'{TfV^{T), 
{AVsYiT) = i'(T)V3(T) = -3i'(T)Vo'(r). 
This yields ( |3.72| ). On the other hand, 
iAVo)^'HT) = 0, 

(ivi)(5)(r) = 5i'(T)y/^^(T) = 5i'(r)(yJ^^-(iyo)^^0(T) = o, 
{Av2)^^\t) = 4:A'{T)V^^\t). 

Since 

V2 = v; - AVi = Fo" - {AVoY - AVi, 
we obtain with (3.66) and ( |5.11| ) that 

v^'^\t) = v^^\t) - {AvoY*\t) - (ivi)(3)(r) = 6i'(r)Vo'(r), 

hence 

{Av2)^^\t) = 24i'(r)Vo'(r). 

On the other hand, 

{AV3)^^\T) = 3A'{T)V^{T) 

= zA{T){Vi{T) + {AV3)'{T)) 

= 3i'(T)(F5(T) + i'(r)F3(r)) 

= mA'{TfVQ{T) (5.18) 

(5.19) 

{AV^YiT) = A'{T)V5{T) = l^A'{T)V^{T). (5.20) 
Gathering together (5.9) and ( 5.14[ )-(5.20), we obtain ( 3.73| ). The proof of Proposition |3.14| is 
complete. 



(5.17) 



where we used (3.51) and (3.71 )-( 3.72). Finally, 



{AVi)"{T) = 2A'{T)Vi{T) = 2A'{T)V^{T) = 30A'(r)Vo(r) 



and 



5.4. Proof of Proposition 3.15, From (3.51)-(3.52 ), we have that 



Thus 



Uo ^ c, Ui = uu - Du^.i - \/i > 1. 

U2{T) = {U[-DUi-Vi){T) 

= {0-{DUoy-V^){T)-V,iT) 
= -D'{T)Uo-2V^{T) 



(5.21) 



where we used successively (5.21), (3.55) and (3.70). 



42 



RODRIGO LECAROSi'2 AND LIONEL ROSIER^ 



Successive applications of (5.21) yield 



Ua{T) 



i=0 



Using (3.66), we obtain that 

3 



i=0 



i=0 



2D'{T){U2{T) + Vi{T)) + D'{T)V^{T) 
-3D'{T){D'{T)Uo + 2yo'(r)) + 2D'iT)V^{T) 
-AD'{T)Vq{T) 



where we used (3.70), (3.74) and the fact that D'{Ty = 0. 



On the other hand, 



i=0 



-2A'{T)V^{T) + 2V3{T) 
-8A'{T)Vq{T) 



by ( [3?7l| ). Combining ([532j)-([534]), we obtain ( [375] ). 



Let us now compute Ue{T). Successive applications of (5.21) yield 

5 

(5-0 , 1/(5-^)1 



i=Q 



We have that 



J2m^)^'-'\T) = J](5 - i)D'{T)ut'\T). 



i=0 



1=0 



Let us estimate the terms uj^'^ ^\t) for i = 0, ...,4. Obviously, Uq'*'{T) = by (5.21), while by 



(4)/ 



(5.22) 



(5.23) 



(5.24) 



(5.25) 



uI^\t) = -{DUo)^'''>{T) - Vq^^'{T) = 0. 



U3) 



(3), 



Next we use (5.21) to obtain successively 



[/^(T) = U4{T) + V3{T), 

U'i{T) = U'^{T) + {DU2)'{T) + Vi{T) 

= UiiT) + V3iT) + D\T)U2iT) + V3{T). 



(5.26) 
(5.27) 
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It follows that 

4 



1=0 



= 3D'{T){U4{T) + 2V3{T) + D'{T)U2{T)) + 2D'{T){Ui{T) + V:^{T)) + D'{T)Ua{T) 
= D'{T){6Ui{T) + 8V3{T)) 

= 2AD' {T){D' {T) + 2A!{T))Vq{T) - 2AD' {T)^ {T)Vq{T) 

= 2AD'{T)A'{T)Vq{T). (5.28) 
On the other hand, using ( |5.15 )-(5.18) and (3.71 )-(3.72 ), we have that 

4 
i=0 

= 6i'(r)Vo'(r) + 2y5(r) + i'(r)y3(r) 

= 33i'(r)Vo(r). (5.29) 



^7J^ follows from ([5^- ([5^. 

Finally, we compute Us{T). We see that 



Us{T) = -Y^[{DU,Y'-^^ + Vl'-'^]{T). 



i=0 



Then 



(7-i) 



(T) 



i=0 



E(7-i)I)'(r)C/,(^-^)(T) 

i=0 

6D'{T)ui^\T) + 5D'{T)U^'^\t) + AD' {T)uf\T) 
+?,D' {T)U'l{T) + 2D'{T)U'^{T) + D'{T)Uq{T). 



Using (3.66), (5.21) and (5.26), we readily see that 

uP{t) = u^^\t) = q. 



Next, successive applications of (5.21) give 



U'l{T) 



uf\T) 



Uq{T) + V^{T), 

U'^{T) + {DUi)'{T) + Vi{T) 

Ue{T) + D'{T)U^{T) + 2V5{T). 

UUT) + {DUs)"{T) + V^'{T) 

{U^{T) + D\T)Ui{T) + 2y5(T)) + 2D'{T){Ui{T) + V^^iT)) 
+V^{T)+A{T)Vi{T). 



(5.30) 
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Thus 



4.D'{T){Uq{T) + 3D'{T)U^{T) + 2V5{T) + 2D' {T)V^{T) + A'{T)V^{T)) 



1=0 



+3L»'(r)(C/6(T) + D\T)Ui{T) + 214, (T)) 
+2D'{T){U^{T) + V,,{T)) + D'{T)Uq{T) 
= D\T)[10U(i{T) + 16V5iT) + 4A'{T)V3iT)] 
= D'{T)[-2AOD\T)A\T)Vo{T) - 330 A' {T)Vq{T) 

+24oi'(r)Vo(r) - i2i'(r)Vo(r)] 
= -io2L>'(r)i'(r)Vo(r). 

It remains to compute Yli=o'^i''' It is easy to see that 

yJ^)(T) = v}^\t) = vi^\T) = 0. 



Successive applications of (3.51) give 

vi;{T) 



vt\T) 



V^{T) + {AV^)'{T) = Vy{T) + i'(r)F5(T), 
Vl!{T) + {AV^)"{T) 
Vj{T) + 3A{T)V^{T), 
V^^\t) + {AV3)(^\T) 
V^{T) + 3A'{T)V^{T) + 3A'{T)V^{T) 

Vj{T) + 6i'(r)y5(r) + 3i'(r)V3(r), 



where we used (5.18). Thus 

7 



^V>^ '\t) = W7{T) + loA {T)V^{T) + 3^ {TfV^iT) 



j=0 



-3M A' {TfV^{T). 



(5.31) 



(5.32) 



Then (3.77) follows from (5.30)-(5.32). The proof of Proposition 3.15 is achieved. 
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